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Abstract 

Charged pions from Ni + Ni reactions at 1.05,1.45 and 1.93 AGeV 
are measured with the FOPI detector. The mean vr^ multiphcities per 
mean number of participants increase with beam energy, in accordance 
with ear her studies of the Ar + KCl and La + La systems. The pion kinetic 
energy spectra have concave shape and are fitted by the superposition of 
two Boltzmann distributions with different temperatures. These apparent 
temperatures depend only weakly on bombarding energy. The pion angular 
distributions show a forward/backward enhancement at all energies, but 
not the = 90'' enhancement which was observed in case of the Au + 
Au system. These features also determine the rapidity distributions which 
are therefore in disagreement with the hypothesis of one thermal source. 
The importance of the Coulomb interaction and of the pion rescattering 
by spectator matter in producing these phenomena is discussed. 

1 Introduction 

This paper presents data on the tt^ emission from Ni + Ni reactions at bom- 
barding energies 1.06, 1.45, 1.93 AGeV. The data were obtained with the FOPI 
detector at GSI and cover almost Air. They are considered a supplement and 
extension of a previous study of the Au + Au reaction at 1.06 AGeV which was 
also performed with the FOPI detector [|l|. The purpose of our new study is to 
systematically investigate the dependence of the pion production on system mass 
and energy in heavy-ion collisions. The present paper is therefore organized in 
close similarity to ref.[Q. 

The pion production rate in heavy-ion collisions increases rapidly with bom- 
barding energy. At 14 AGeV the number of pions measured in central collisions 
has reached the number of baryons contained in the participant 0. At energies 
where the number of mesons is still smaller than the number of baryons, calcula- 
tions based on the BUU model [Q have shown that the number of pions reaches 
its final value, which is the one measured, only at very late stages of the reaction. 
At earlier times pions are partly absorbed by the nucleons with the result that 
a considerable part of the baryons consists of nucleon resonances, primarily of 
the A(1232) resonance. For the observed properties of pions the concept of their 
freeze-out time and the complex evolution of the baryon distribution are impor- 
tant considerations. Pions are emitted during the total time of the reaction and 
thus their phase space distribution bears the signature of the complete reaction 
history. The early measurements of the pion energy spectra ||^ showed that 
these are quite different from the corresponding proton spectra. Compared to 
pure Boltzmann spectra with fixed temperature T the former are, at around 1 
AGeV bombarding energy, concave shaped with a mean temperature of about 
T = 70 MeV, whereas the latter are convex shaped with a mean temperature 
of about T = 100 MeV. This difference was first interpreted as due to the 
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different source sizes at freeze-out time, later it was realized [0 that pions mainly 
originate from the decay of the A(1232) resonance, and that their spectral shape 
is strongly influenced by the decay kinematics. Therefore pions are very good 
indicators for the excitation of such resonances. On the other hand, besides the 
decay channel A ^ + nn, also the competing channel A + ^ + de- 
termines the fate of the resonances in nuclear matter. The time evolution of the 
nuclear matter distribution is, in addition to the decay kinematics, the second 
decisive factor for the pion phase space distributions. 

The question of how these processes change with mass and energy of the 
colliding nuclei motivated the present investigation. The comparison with other 
reactions particularly include data from ref. and from the work of Harris et 
al.p]. Recently corroborative data from other experiments on pion production. 



performed by the TAPS H and KaoS |10| collaborations at GSI have become 



available. Besides these published results, the more detailed GSI reports M IjlU 



and conference proceedings, as e.g. ref. |TT|, are of particular interest for the 
present investigation. 

The paper is organized as follows: The section 2 gives a short summary of the 
experimental procedures. This part is shortened to the absolute necessary since 
a more detailed account of these procedures was published in ref. ^ . The section 
3 presents the experimental results, frequently data from the Au -|- Au reaction 
are included to illuminate the dependence on system mass. The summary and 
discussion of the results on pion production obtained so far with the FOPI 
detector are contained in section 4. 



2 Experimental procedures 

The reaction Ni + Ni was studied at nominal beam energies of 1.06, 1.45 and 
1.93 AGeV. The ^^Ni beam was accelerated by the UNILAC / SIS accelerator 
combination of the GSI /Darmstadt. The duty cycle was 75% with a spill length 
of 4 s. The average beam intensity corresponded to 5 ■ 10^ particles per spill. 
The target consisted of a ^^Ni foil of 270/im thickness. This corresponds to an 
interaction probability of 0.5%. The average energy loss of the ^^Ni beam in this 
foil amounts to 0.005 AGeV and was neglected. 

The particles produced by the Ni -|- Ni reactions were detected by the FOPI 
detector which is a modular detection system with almost 47r coverage. The 



FOPI detector is described in [|T^ |jT3[ . Of particular importance for the present 
investigation are the central drift chamber CDC which is mounted inside the 
superconducting magnet with a solenoidal field of 0.6 T strength, and the for- 
ward plastic scintillation wall PLA. Both detector components have complete 
cylindrical symmetry, the CDC covers the polar angles from 32*^ to 150*^ in the 
laboratory frame and allows pions to be distinguished from baryons. The parti- 
cle multiplicity measured by the CDC is labelled ncoc, when pions are excluded 
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we use as quantity the measured baryonic charge ^cdc- '^^^ PLA covers the 
polar range from 7^ to 30°, it allows the separation of particles only accord- 
ing to their charge number Z. The measured particle multiplicity in this case 
is labelled tt-pla- Similarly the total particle multiplicity is labelled nxoT with 
^TOT = ncDC + ^PLA- Because of the inability of the FOP I detector at the time 
of the experiment to identify pions for polar angles smaller than 32° the solid 
angle of the pion identification amounts to around 8.2 sr. But pions emitted into 
the uncovered area of the complete solid angle can be reconstructed by employing 
the symmetry relation 

/(e,<l>) = /(7r-0,7r + <l>), (1) 

where / may be any angle dependent observable and B, $ = refer to the cm 
frame. This relation is valid only for symmetric reactions. Including the remain- 
ing detection and geometrical thresholds it is estimated that after symmetrization 
around 90% of the pion momentum space becomes accessible. 

In the present study we shall use the identical conventions which were used in 
ref.||l|. All quantities which refer to the laboratory(Za6) or target frames will be 
labelled by small letters. The center of mass(cm) or fireball frames are generally 
characterized by capital letters. Exceptions are those quantities, like the trans- 
verse momentum, which remain unchanged under the transformation of frames. 
The magnitude of the transverse momentum pt is given in units of pt°'' where the 
index *^°) indicates a normalization by the factor {A ■ -Pproj/^proj)"^ and where A 
is the particle's mass number. Similarly the rapidity Y in the cm system is nor- 
malized by the factor (l^proj)"^ and this quantity is labelled Y^^\ In the present 
experiments the normalization factors for the pion transverse momentum have 
the values 9.47, 8.09, 7.02(GeV/c)-i for 1.06, 1.45, 1.93 AGeV bombarding energy, 
equivalently the normalization factors of the rapidities are 1.44, 1.26, 1.12. 

The data from the drift chamber CDC were analyzed by using three different 
tracking algorithms as described in ref. [Q. The purpose of this threefold analysis 
is predominantly to obtain a measure of the systematic errors involved in the track 
finding and particle identification. Because of the smaller particle multiplicities, 
compared to the Au + Au reaction, the systematic errors are reduced in the 
present case. They are largest for the pion multiplicities where, for large impact 
parameters, they become as large as 10%. For other quantities, like the pion 
mean kinetic energies or the widths of their rapidity distributions, they always 
remain smaller than 5%. In [|l| the uncertainty of measuring the pion transverse 
momenta from the curvature of the pion tracks was determined to increase form 
4% to 7% in the range < pt < 1000 MeV/c. From the accuracy of ±2° with 
which the forward CDC limit i} = 32° is identified in the particle spectra one 
obtains a similar systematic uncertainty of 7% with which the pion lab momenta 
can be measured in forward direction. In this paper we generally quote only one 
error which is the larger one of the statistical and systematic errors. Normally 
the latter determines the error. 
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Table 1: Lower multiplicity boundaries nTOT{TM5), cross sections crroac, and 
deduced radii tq for the Au + Au (first row) and the Ni + Ni reactions 



energy 


nTOT(TM5) 




^0 


(AGeV) 


lower limit 


(barns) 


(/m) 


1.06 


146 


5.50 ± 0.50 


1.14 ± 0.05 


1.06 


60 


2.19 ± 0.24 


1.08 ± 0.06 


1.45 


64 


2.44 ± 0.29 


1.14 ± 0.07 


1.93 


68 


2.35 ± 0.21 


1.12 ± 0.05 



3 Experimental results 

3.1 Charged particle multiplicities 

The multiplicity distributions of charged particles measured with the forward 
scintillation detector PL A and the central drift chamber CDC vary only weakly 
for the three chosen bombarding energies. These spectra rise for small multiplic- 
ities riTOT, their plateaus for large multiplicities are less well developed than in 
case of the heavier Au + Au reaction, c.f. Fig. 6 of ref.[|l|. It is common practice 
to subdivide the multiplicity distributions into 5 bins TMl to TM5 where the 
upper bin TM5 has its lower multiplicity boundary at half the plateau value, and 
where the other 4 bins have equal multiplicity widths. These bins select event 
classes of certain impact parameter ranges, with TMl corresponding to the most 
peripheral reactions, and TM5 to the most central reactions. The values which 
define the lower multiplicity boundaries of the TM5 bins are listed in table 1. The 
slight increase of this value with bombarding energy is mainly due to an increase 
of the charged pion multiplicity. The 'minimum bias' situation is defined by the 
condition that events are not selected according to their total charged particle 
multiplicity nxoT- 

After the integration of the multiplicity distributions over their total range, 
and after corrections because of acceptance losses for multiplicities rixoT < 8, 
the total cross sections a^ea.c for the 3 bombarding energies are deduced and 
listed in table 1. This table also includes the corresponding information from 
the Au + Au reaction. These cross sections define the nuclear radii tq via 



ro = 0.5y4~^/^y (Treac/TT whcre A is the mass number of the target respectively pro- 
jectile, and where it is assumed that both nuclei have sharp surfaces. The values 
of ro are listed in table 1. They are for all systems studied close to ro = 1.12/m 
and therefore slightly smaller than the values deduced from electron scattering 
It was argued in ref.[jl| that the agreement with electron scattering is im- 



proved when the smooth nuclear surface is taken into account. 

It was already pointed out that the multiplicity nxoT is related to the impact 
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50 100 50 100 50 100 



part part part 

Figure 1: The dependence of the tt^ multiphcities on the number of nucleons 
Apart in the fireball. The full and dotted lines present the extrapolations from 
the results of Harris et al.(see text) 



parameter b and thus to the size of the nuclear fireball, expressed by the number 
Apart of participants. To deduce the quantitative relation between titot and Apart 
requires to parametrize the nuclear density distribution by a Woods - Saxon form 
factor as was outlined in ref. Applying the identical procedures and also using 
the identical global quantities nTOT, ''^pla, -^^cdc number of charged pions n^ri 
can be obtained as function of Apart- We show the dependence of n^ri on Apart 
in Fig.l, where the measured pion multiplicities were used to calculate the mean 
multiplicities at regular intervals of Apart- The measured pion multiplicities were 
in all cases multiplied with a factor 1.1 to correct for the remaining detection 
inefficiencies [Q. This 10% correction of the pion multiplicities is predominantly 
caused by the = 32° cut for pion transverse momenta pt < 0.6, c.f. Fig.3. Its 
size was estimated using the IQMD model [1^ or a thermal model in which pions 
are isotropically emitted with the measured mean kinetic energies. Owing to these 
models which both do not describe the experimental observations completely 
0, the accuracy of this correction cannot be expected to be better than 50%. 
Within this limit the correction does not change with system mass or energy or 
impact parameter, since the pion momentum distributions do not change very 
much, as discussed later. The horizontal errors present the estimated accuracies 
with which Apart can be obtained from titot, 't^pla, ^cdc' vertical errors are 
predominantly systematic and present the uncertainties which were estimated 
by using different trackers. In addition the Fig.l displays the expected pion 
multiplicities which one obtaines by extrapolating the Fig. 2 of |Q] into the Ni + Ni 
system. The comparison with the expectation demonstrates that the dependence 
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Table 2: Expansion coefficients a\ for the expansion of (Apart) in powers of 
Apart for the Ni + Ni reactions 



energy 


TT 


71+ 


(AGeV) 


af) ■ 10^ 


af) ■ 10^ 


af ) • 10^ 




1.06 


4.22 ±0.89 


0.41 ±0.71 


5.76 ±0.76 


-1.96 ±0.67 


1.45 


7.30 ± 1.26 


-0.41 ±0.94 


8.73 ± 1.22 


-3.02 ±0.95 


1.93 


11.4 ± 1.53 


-1.41 ± 1.18 


12.5 ± 1.21 


-4.12 ± 1.06 



of r?,,r± on Apart IS not hnear in Apart, but that the mean charged pion multiphcities 
should be close to the results obtained by Harris et al. P] when extrapolated into 
the Ni + Ni system. In order to account for the measured non-linearities n^i was 
described by a power series of second order in Apart, i-e. 

^7r(Apart) = (^1 ^Apart ± 0^2 ^^part' (2) 

where the index ^"^^ indicates the order of the expansion. The expansion coeffi- 

(2,) 

cients a- are listed in table 2 for the three energies used. In a way described 
in ref. the average number of pions per average number of participants can be 
deduced from the ap-* via 



"1 



(2) 

a? + ^Ao, (3) 



< Apart > 2 

where Aq is the total number of nucleons. The results are shown in table 3 which 
includes also the Au + Au reaction. Notice that these values depend on the order 
of the expansion. If only a 1. order expansion is chosen, as was additionally done 
in 0, the vr" multiplicities in the present case were not effected, but the 7r+ 
multiplicities would be reduced by 15%. This would yield a similar increase of 
the 71" to 71+ ratios, but an 8% reduction of the total pion multiplicities. The 
total pion multiplicity < 71-^ > is calculated using the isobar models prediction 

m 



{5N^ + NZ) ■ {N^ + mZ ± Z^) : (5^^ + NZ) (4) 

where in this case iV, Z are the numbers of neutrons, respectively protons in the 
target or projectile. The eq.(4) holds if the pions are emitted by a resonance 
with isospin / = 3/2, and yields for the n^- to ?t,^+ ratio a value i?^*-* = 1.12. 
This is in very good agreement with the experimental mean values of the ratio 
R'^'^ =< n^^ > / < > also shown in table 3, whereas in case of the Au ± 
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Harris extrapolation J 
• Au + Au .-'^ 
■ Ni + Ni i.-* J 

♦ 

/ ♦ ] 

I I I L>'* I I I I I I I I I I I I I I I I 
0.25 0.5 0.75 1 " 0.' 0.7 0.3 0.4 0.5 

normalized Ap^^ cm bombarding energy (GeV) 

Figure 2: Left: The dependence of the normahzed 7r~ to 7r+ ratio R^yR^^ on 
the scaled number of participants 74pait/^o foi' different reactions and bombarding 
energies. The dotted hue corresponds to R^^ = R^^\ Right: The comparison be- 
tween the measured mean number of pions per participant and the extrapolation 
from the results of Harris ct al.(scc text) 

Au reaction one finds a considerable disagreement between the model prediction 
R^^ = 1.95 and the experimental value. If resonances with isospin 1 — 1/2 were 
also involved the ratios in eq.(4) would be replaced by 

n^- : n^o : = 2Z : {N + Z) : 2N , (5) 

which is smaller than the experimentally observed value. One is, however, not 
allowed to conclude that this proves a contribution from 1 = 1/2 resonances in 
case of the Au + Au reaction, since eqs.(4,5) are vahd only for first-generation 
decays and do not include multiple generations of resonances. 

This result appears rather puzzling because the non-linear behavior of the 
pion multiplicities n.n-± nevertheless yields a functional dependence of the ratio 
R^^^ on Apart which is very similar for all reactions studied. In order to demon- 
strate this similarity we have normalized R!^^ to R^^ and Apart to Aq, and these 
normalized quantities are plotted in Fig.2. After the normalization the 7r~ to tt"*" 
ratio increases almost linearly with the normalized Apart, for small values of Apart, 
i.e. large impact parameters, one finds an apparent deficit of vr" whereas for cen- 
tral collisions there appears to be a deficit of tt"*". Notice that the dependence 
of R^^ on Apart does not disappear when Apart is determined by means of the 
charged particle multiplicity but employing a different method, since a different 
method will not change the values of R^^^ although it might change their exact 
dependence on Apart- Nevertheless, the condition of total charge conservation 
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Table 3: Average number of pions 71^^ per average number of participants Apart 
and average n^^- to n.,^+ ratio for Au + Au (first row) and Ni + Ni reactions. 
Only systematic errors are shown, the statistical errors are in the order of 10% 
of the systematic errors 



energy(AGeV) 


<^~>(%) 

< Apart > ^ ' 




< Apart > ^ ' 




1.06 


3.08 ±0.37 


1.82 ±0.43 


7.26 ±0.84 


1.69 ±0.26 


1.06 


4.46 ±0.51 


4.62 ±0.41 


13.62 ± 1.39 


1.19 ±0.07 


1.45 


7.06 ±0.80 


6.98 ±0.75 


21.05 ±3.29 


1.04 ±0.17 


1.93 


10.55 ±1.01 


10.07 ±0.79 


30.92 ±3.83 


1.05 ±0.13 



can induce a correlation between the measured particle multiplicity and the ratio 
R^) which will be discussed in more detail in section 4. 

In Fig. 2 we have also plotted the dependence of the average number of all pions 
per participant on the bombarding energy. As expected the Ni ± Ni system 
closely follows the systematics obtained by Harris et al.^, whereas the pions 
from the Au ± Au system are reduced in number by almost a factor 2. It was 
pointed out in ref . [T[] that this reduction is confirmed in measurements of the vr^ 
production by the KaoS collaboration and in case of the 7r° by the TAPS 
collaboration |[ll|. Both collaborations use detectors which have a solid angle 



smaller than 47r and therefore require the extrapolation from the measured range 
into complete solid angle. This extrapolation is not simple since the angular 
distributions of pions are non-isotropic, c.f. section 3.2, and since the anisotropy 
depends on system mass. But the observed reduction is of a size that even 
a better extrapolation procedure will not reconcile these results from different 
experiments with the expectations one obtains by extrapolating from the light- 
mass to heavy-mass systems. 



3.2 Phase space distributions of charged pions 

In the previous investigation of the Au ± Au reaction the phase space distribu- 
tions of vr^ were found to be non-thermal, requiring at least two temperatures 
and a non-isotropic angular distribution for their parametrizations. Similar ob- 
servations are made in the present case of the Ni ± Ni reaction as will be shown 
in detail in the following subsections. 

The Fig. 3 shows a total view onto the vr" and vr"*" momentum space distri- 
butions under minimum bias condition in case of the Ni ± Ni reaction at 1.06 
AGeV bombarding energy. Also shown are the curves of constant momentum 
P = 0.2, 0.4, 0.6 GeV/c to illustrate the deviations from a pure thermal behavior, 
and as broken curve the i} = 32° geometrical limit of the CDC. Pions at smaller 
angles were not measured but obtained by employing the reflection symmetry 
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cm rapidity ' cm rapidity Y 



Figure 3: The momentum space distribution under minimum bias condition for 
7r^ and 7r+ at 1.06 AGeV bombarding energy. The full curves are displayed for 
pion momenta P = 0.2,0.4,0.6 GeV/c, the dotted curves display the i} = 32° 
forward boundary of the CDC. The invariant cross section (Pa/{ptdptdY^^^ is 
displayed in arbitrary units on a logarithmic scale 

around Y^'^^ = 0. The Fig. 3 should be compared to the equivalent Fig. 10 in 
ref.||l[| to see that differences between the pion emission from the Ni + Ni and 
Au + Au are only small. Nevertheless such differences exist, and their detailed 
exposure is the subject of the following subsections. 

3.2.1 Kinetic energies 

The kinetic energy distributions of tt~ and vr^ from the Ni + Ni reactions have, 
at all bombarding energies studied, concave shape, and therefore two Boltzmann 
distributions with different temperatures, a low temperature T\^t^ and a high tem- 
perature Th,7r, are needed to fit these distributions in the frame of a thermal 
picture. This is similar to the Au + Au reaction presented in ref.|jl|], where 
the Fig. 12 showed a representative example of the measured distributions at 
B = 130°. Similar to the Au + Au reaction is also that the extracted values for 
Th,7r appear to increase for emission angles G > 140°. The angular dependence 
of the low and high-temperature values is plotted in Fig. 4. The values show 
considerable fluctuations, the fluctuations are inherent to our procedure which 
flts a sum of two exponentials to a distribution which only slightly deviates from 
a pure exponential behavior. The temperature increase at backward angles is 
seen most strongly for the tt"*", in order to illustrate the size of the increase the 
Fig. 5 compares the tt~ spectrum at G = 90° to the vr"*" spectrum at G = 150°, 
where the pions were measured at 1.06 AGeV bombarding energy. The Fig. 4 also 
reveals that Ti 7r+ is systematically larger than Ti^n-, whereas for G < 140° the 
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Figure 4: The apparent pion temperatures Ti and Th for different cm angles 
G. The hues display the the mean values in the range 85° < < 135° 

Table 4: Apparent pion temperatures Ti and Th^vr and mean ratio of the low- 
temperature to high-temperature component -Rt for Au + Au (first row) and Ni 
-|- Ni reactions. The error shown is the larger of the statistical and systematic 
errors for a given entry 



energy (AGeV) 


Ti,.-(MeV) 


Ti,.+ (MeV) 


Th,.(MeV) 




1.06 


42.2 ±2.7 


49.4 ±3.7 


96.4 ±5.1 


0.75 ±0.36 


1.06 


47.0 ± 1.9 


52.5 ±3.6 


93.5 ±8.0 


1.30 ±0.81 


1.45 


51.7± 1.9 


54.3 ±3.3 


103.8 ±4.9 


1.15 ±0.37 


1.93 


53.1 ±5.7 


56.6 ±4.4 


111.7±9.5 


1.06 ±0.80 



fits to the energy distributions are consistent with the assumption that T}^^^+ is 
equal to Th^^r- . The temperature values obtained from such fits depend on the fit 
range. In our analysis the vr" kinetic energy spectra were fitted in the range from 
50 to 750 MeV, whereas the fit range for the vr"*" kinetic energy spectra varies 
with angle 0: for 6 = 90° the range is from 50 to 375 MeV, the upper limit 
has increased to 750 MeV for = 130°. As pointed out in ref.|l| this angular 
variation of the fit range is necessary to separate vr"*" from protons in the analysis. 

In accordance with the analysis of the pion energy spectra performed in ref. |ll], 
and because of the increase of the angular cross section da/dQ for > 140° (see 
below), we have extracted the average pion temperatures only from pions emitted 
into the angular range 40° < < 140°, the values are listed in table 4 and shown 
in Fig.4. The values of Th^^r agree with an independent analysis of the present 



11 



5- ^""^ 

LU 

■o 



■o 

* 

LU 10' 
D. 

1 









• 7i' 0=90° 


•a 

^ \ 

r ^ 













0.2 0.4 0.6 0.8 

cm kinetic energy (GeV) 



Figure 5: The kinetic energy spectra of vr and tt^ observed at 1.06 AGeV bom- 
barding energy at different cm angles G 



data based on the transverse momentum spectra [lYj, but they are significantly 



larger 20) MeV than those obtained by the KaoS collaboration [jlO[. In the 
latter case the 7r+ momentum spectra, measured at ^ = 44° ± 4P in the lab frame, 
were analyzed. Whether this difference in the method or different fit ranges are 
responsible for the discrepancy is not clear. The table 4 also contains the ratio 
Rt, which is the ratio between the low-temperature component Ji^tt and the high- 
temperature component /h,7r, and it contains the equivalent results from the Au + 
Au reaction. The intensities Ji^^r and /h,7r were calculated by integration over the 
Boltzmann distributions with temperature Ti ^r, respectively Ty^^n- The differences 
between Ti .n.+ and Tj .^.-, which amount in the average to 3.9 ± 1.5 MeV in case 
of Ni -|- Ni, but to twice this value in case of Au + Au, has its origin in the 
Coulomb repulsion, respectively attraction between the pions and the positively 
charged baryonic matter. Therefore the ratio is influenced by the Coulomb 
interaction, and to compensate this effect the values of R^ quoted in table 4 are 
the mean values between vr" and 7r+. Notice that the large errors assigned to the 
values of Rt are mainly systematic and due to the fact that Rt changes more 
with angle for O < 140° than the temperatures. 

Both pion temperatures, the low Ti temperature and the high T^^^^ tempera- 
ture, increase slightly with beam energy. Similarly also the relative contribution 
of pions with high temperatures increases, but this component never exceeds the 
contribution of low-temperature pions in the cases of the Ni + Ni reactions. This 
is different in the Au + Au reaction where already at 1.06 MeV bombarding 
energy the high-temperature component is stronger than the low-temperature 
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Figure 6: The mean n.^— to 71^^+ ratio in the range 85'^ < < 135° as function of 
the cm kinetic energy. The dotted hne corresponds to equal numbers of 7r~ and 



TT ' 



component. 

The use of the term 'temperature' in discussing the kinetic energy distri- 
butions of charged pions should not be confused with the true thermodynamic 
temperature describing unordered motion. Even in the case that pions thermalize 
in nuclear matter their finite mean free path adds a zero point energy to their 
thermal energies because of their small mass. The size of this quantum effect 
is in the order of 25 MeV, its general dependence on mass and localization of 
the pions was calculated in ref.|]l8|. In our present analysis the temperature is 
solely used as a fit parameter to reproduce the concave shapes of the energy 
distributions. Its small sensitivity to the bombarding energy indicates, that the 
freeze-out conditions of pions do not depend strongly on the initial conditions, 
on the other hand they are not completely independent of these conditions. It 
it obvious that to make the differences visible, the experimental accuracies have 
to be accordingly high, and the fit conditions have to be defined in a consistent 
way. 

For the sake of completeness and in accordance with the deduced pion tem- 
peratures we display in Fig. 6 the ratio of the kinetic energy spectra of vr" to vr"*" 
under minimum bias conditions, where only the cm angles 90" < 6 < 135° are 
included. For the three bombarding energies studied the dependence of this ratio 
on the pion kinetic energy changes only very little and is, except for the smaller 
values at small kinetic energies, similar to the Au + Au case, c.f. Fig. 13 in ref. |l[]. 
Notice in particular that the intersection with one occurs in all systems at kinetic 
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Figure 7: The pion angular distributions for peripheral (open symbols) and cen- 
tral (full symbols) reactions 



energies less than 300 MeV. The global dominance of the vr" over the 7r+, ex- 
pressed by R^^ > 1 is therefore caused by the more abundant low-energy pions, 
at larger energies this relation appears to be reversed. The observation that the 
7T~ yield becomes smaller than the vr^ yield for kinetic energies above 100 MeV 
is linked to the reactions with large impact parameter, c.f. Fig.l and Fig. 2. In 
the case of central (TM5) collisions the tt^ to tt'^ ratio approaches the value one 
for kinetic energies above 200 MeV. Accordingly the shape of the kinetic energy 
spectra of pions, i.e. the apparent pion temperatures, do not change significantly 
with impact parameter, but the yields do. 



3.2.2 Angular distributions 

The angular distributions of charged pions were found in previous experiments 
to be non-isotropic with a pronounced enhancement at forward and backward 



angles [0. This feature is similarly observed in the present Ni + Ni reactions 
at all energies, c.f. Fig. 7. In the Fig. 7 the measured angular distributions da/dQ 
are displayed under the selection criteria of peripheral (TMl) and central (TM5) 
collisions. All measured pions with kinetic energies above 50 MeV are included. 
Notice that the reversal in the sizes of the vr" and vr"*" cross sections between 
TMl and TM5 is a phenomenon similarly seen in Fig.l or Fig.2. Since in the 
elementary N + N^N + N + tt reaction pions are also preferentially emitted 
into forward and backward angles, one possible explanation for the anisotropy 
found in nucleus - nucleus reactions is based on first-chance collisions in the 
dilute surfaces of the colliding nuclei. Alternatively the rescattering of pions in 
spectator matter yields a similar angular anisotropy. It was argued in lel.M 
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that this latter mechanism could also explain the increase of the apparent pion 
temperatures at forward and backward angles. Contributions from spectators 
are always assumed to be present in the phase space distributions of baryons, 
and a similar anisotropy in the angular distributions of protons under minimum 
bias conditions is readily interpreted as due to spectator decay. Indeed, the 
forward/backward enhancement changes more strongly with impact parameter in 
the Au + Au reaction than in the Ni + Ni reactions, and this change is different 
for 7r~ than for vr"*". In the first reaction the anisotropy o"(150'^)/cr(90°) increases 
between TM5 and TMl by a factor 1.9 ±0.2 for vr" and 1.4 ±0.2 for vr"*", whereas 
in the second reactions the corresponding factors are 1.3 ± 0.2 and 1.0 ± 0.2 
independent of the bombarding energy. The differences in the surface to volume 
ratios for Au and Ni, and the effects of the Coulomb attraction, respectively 
repulsion may be responsible for this behaviour of the pion angular distributions. 
It suggests that even in central collisions a large part of the nucleons is not 
stopped in Ni ± Ni reactions. 

The pion angular distributions reveal another peculiarity at angles around the 
transverse direction = 90*^. Whereas pions from central collisions of the Au 
± Au system display a clearly discernible enhancement of the emission into this 
angular range, c.f. Fig. 17 of ref. |l[], this enhancement is not observed for pions 
from the central Ni + Ni collisions at all energies studied. Earlier measurements 
of the pion angular distributions from light-mass systems gave a similar result 
0, i.e. a flat da/dQ at angles around O = 90°. Indeed, the transverse pion en- 
hancement in the Au ± Au central collisions constitutes the most direct evidence 
that the pion phase space distributions change with system mass at bombarding 
energies around 1 AGeV. The cause for this change is not understood. It is remi- 
niscent of the predicted hydrodynamical flow in heavy systems when the nuclear 
compression approaches values of three times normal nuclear density [|19|. But 
it is not clear why this is observed in the pion channel when it is most likely not 



present in the baryon channels |2^. With this respect it should be remembered 
that apparent pion flow phenomena are quite often interpreted as caused by pion 
absorption in the nuclear environment ||2l| [E^ . 



3.2.3 Rapidities 

The forward/backward anisotropy seen in the pion angular and temperature dis- 
tributions also determines the shape of the pion rapidity distributions. These 
are displayed in Fig.8, they are obtained from the complete phase space dis- 
tributions shown in Fig.3 for a special case, by a projection onto the rapidity 
axis with pt thresholds of pt > 85 MeV/c (1.06), 100 MeV/c (1.45), 115 MeV/c 
(1.93), where the number in brackets specifies the bombarding energy in units 
of AGeV for which these pt cuts were applied. The dotted curves in Fig.8 show 
the expected rapidity spectra, if pions were isotropically emitted with the two 
temperatures Ti ,r and T^^t, determined in subsection 3.2.1. It is obvious that the 
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Table 5: Parameters of the three-source fits to the pion rapidity spectra under 
minimum bias condition for the Au + Au (first row) and Ni + Ni reactions. Notice 
that the values for the widths a are presented in units of Y^^'^ . The numbers in 
parenthesis give the systematic error in the last 2 digits. The statistical error is 
of the order (05) for -Ry.tt and (02) for cro,7r, o^i,-k 



energy(AGeV) 




O"0,7r- 


0"l,7r- 




O"0,7r+ 


0"l,7r+ 


1.06 


1.68(21) 


0.45(04) 


0.40(03) 


1.52(14) 


0.45(03) 


0.42(03) 


1.06 


1.21(05) 


0.48(01) 


0.45(01) 


1.16(07) 


0.46(01) 


0.47(02) 


1.45 


1.54(15) 


0.49(01) 


0.42(02) 


1.30(10) 


0.48(01) 


0.44(02) 


1.93 


2.68(13) 


0.51(01) 


0.38(03) 


2.43(34) 


0.53(02) 


0.39(03) 



rapidity distributions deviate from these expectations, the deviations become less 
visible with increasing energy. In ref.||l[ it was assumed that the excess intensity 
at rapidities around \Y^'^'>\ = 1 is due to rescattering by spectator matter, and 
the strength of this process was estimated by fitting the measured spectra with 
three Gaussian distributions located at = — 1,0,+1, where the Gaussians 
at Y^^^ = ±1 have to be identical for a symmetric reaction. This procedure 
was repeated in case of the Ni + Ni reaction. As in the Au + Au case the 
fit parameters, i.e. the width (Ty.tt of the Gaussians in units of Y^^\ and their 
relative intensities Ry,tt = Io,Tr/ {h,Tr + I-i,n) do not change by more than 20% 
with impact parameter. Therefore the table 5 presents the values of ay^n and 
i?Y,7r only under minimum bias condition, in Fig. 9 the ratio Ry^n is plotted for 
the different impact parameters and for the three bombarding energies of the 
Ni + Ni reaction. In addition we show in Fig.9 the widths aY,iT as functions of 
bombarding energy this time in units of the unnormalized rapidity Y. The ratio 
RY,Tr and the width ao^jr increase with bombarding energy, the widths <Ji^^,a_i^T, 
stay constant within their errors. This indicates that the number of pions being 
rescattered decreases with energy, and that the separation into 3 pion sources be- 
comes increasingly ambiguous. At energies above 2 AGeV the rescattered pions 
from the spectator sources are almost completely masked by the non-rescattered 
pions from the central source. Note, however, that the pion angular distributions 
still show clear evidence for the presence of the rescattering mechanism at even 
the highest bombarding energy. Compared to the Au + Au reaction at 1.06 
AGeV the Ni -|- Ni reaction at the same energy has a stronger rescattering signal, 
and this rescattering contribution to the pion rapidity spectra varies less with 
impact parameter than observed in the former reaction. These observations are 
the prominent indicators from the pions momentum space distributions that the 
nuclear corona phenomena and the related formation of nuclear spectators are 
different in the Ni + Ni and Au + Au systems. And it might provide the basis on 
which one might eventually find the cause for the transverse pion enhancement 
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Figure 8: The pion rapidity spectra at different bombarding energies. Tfie fuli 
curves sliow tlie fit by 3 Gaussian distributions (see text), tlie dotted curves are 
tfie tfiermal expectations assuming two temperatures Ti ,r and Th^^r 



in central Au + Au collisions. 



4 Summary and discussion 

The results of this investigation should be discussed in relation to the study 
of pion production in the Au + Au system which contains almost four times as 
many nucleons. The total number of pions produced at a given energy is therefore 
larger in the heavier system, but when normalized to the number of participants 
the Ni + Ni system produces almost two times more pions at 1.06 AGeV than 
the Au + Au system. Furthermore the Ni result is in agreement with equivalent 
results from the Ar + KCl and the La + La systems studied by Harris et al.]^. 
The surface to volume ratio varies from the Ni to the La system (1.34) more 
than from the La to Au system (1.12), and thus it is difficult to attribute the 
reduction of the vr^ yield in the Au + Au case to an enhanced pion absorption 
in the larger fireball volume. To decide whether or not absorption is responsible 
for the reduced tt^ yields, the independence of the normalized number of pions 
for system masses below Aq = 280 should be confirmed, e.g. for the Aq = 200 
system. Considering the fact that the pion angular distributions depend on the 
system mass, it is not sufficient to do this test at midrapidity but it should be 
conducted over the complete solid angle. 
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Figure 9: Left: The impact parameter dependence of the intensity ratios Ry 
between pions from the fireball and pions from the spectators. Right: The energy 
dependence of the widths uy of the pion rapidity distributions from the fireball 
and the spectators under minimum bias conditions 



The alternative explanation for the reduced production of pions in the Au + 
Au system might be an enhanced flow of energy into collective motion which was 



recently observed in central collisions, although at lower energies |20]. This will 
decrease the energy available to excite baryon resonances. At energies around 1 
AGeV pions are considered to mainly originate from the decay of heavy baryon 
resonances, where the A(1232) resonance dominates the pion production. Nev- 
ertheless, theoretical studies have shown that resonances up to a mass of 1.9 
GeV/c^ have to be included to explain the pion spectra measured in Au + Au 



reactions |]23|. The excitation of such resonances is seen in the concave shape of 
the TT kinetic energy spectra. The apparent higher temperature Th^^r is a signal 
for the presence of resonances above the lowest lying A (1232) resonance, as will 
be discussed in more detail in a forthcoming paper [|^. With this respect it is 
interesting to notice that Th^jr is larger in the Au + Au reaction than in the Ni 
+ Ni reaction at the same energy, and that also the contribution of the high- 
temperature component is larger in the Au than in all the Ni cases studied. This 
is a puzzling result, since an enhanced excitation of the higher baryon resonances 
should lead to an increase in the number of pions, because the 27r decay channel 
becomes more important. The conflict can be resolved in a simple way, if baryon 
resonances are only excited via N + N collisions in the early compression phase 
of the collision. In the later expansion phase their decay by pion emission will 
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produce second and subsequent generations of resonances, not only in the fireball 
but also in the cold spectator matter. It is reasonable to assume that the number 
of resonance generations is larger in the heavier system, and that the competing 
N + B N + N channel causes a loss of pions. In this case the total pion 
reduction would be proportional to a power law with the number of generations 
in the exponent. 

Coulomb effects are present in the pion spectra, most indicative is the differ- 
ence between the values of the low temperature Ti ^^- and Ti^+. This difference 
is only half as large in the Ni + Ni reaction than in the Au + Au reaction, and 
appears to be independent of bombarding energy. On the other hand the vr" to 
TT"*" ratio also depends on the size of the fireball, expressed by the number of par- 
ticipants Apart- And this dependence, when properly normalized, is found to be 
independent of the system mass or energy. It has a universal dependence with too 
few 7r~ at peripheral reactions, and too few tt"*" at central reactions. The latter 
might be due to the condition that the total charge is conserved. This condition 
and the requirement that for minimum impact parameter the measured particle 
multiplicity attains its maximum value selects events with large numbers of vr", 
i.e. with large R^^^ values. It is easy to see that the two conditions 

Z^""' + {1 + R^^^)n^+ = max 
Z^""' + (1 - R^f)n^+ = const 

are the better fulfilled the larger R^\ Of course, R^^^ in these equations cannot 
grow unlimited, it is bound by isospin and baryon conservation. But it is conceiv- 
able that it becomes largest when the number of measured particles is largest. In 
the case of peripheral reactions the pion rescattering process in spectator matter 
is probably non negligible in determining the vr" to 7r+ ratio. The impact param- 
eter dependence of the pion angular distributions is suggestive: In both systems, 
i.e. Au -|- Au and Ni -|- Ni, the forward/backward enhancement in peripheral 
reactions is stronger for vr" than for tt"*". Interpreted as due to rescattering, the 
rescattering process in spectator matter would then be stronger for tt" than Tr"*", 
and the Coulomb attraction between the tt" and the positively charged specta- 
tors offers a possible explanation for this preference. The preferred absorption of 
vr" in spectator matter, without reemission because of the N + B N + N chan- 
nel, would also mean that pion absorption is an important process to heat the 
spectators. In addition the preferred absorption of negative pions should cause 
the decay products of the spectator to become more neutron-rich. The effects 
the Coulomb potential has on the ratio R^^ was recently studied in a theoretical 
paper by Teis et al. [|^. Whereas its dependence on the pion kinetic energies, 
c.f. Fig. 6, is quite well reproduced, the dependence on the number of participants 
Apart deduced from the measured particle multiplicities is not studied in similar 
detail. 

With respect to the pion 'temperatures' Ti^^ and Th^Tr listed in table 4, we 
want to reiterate that these should be only considered as fit parameters to the 
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concave shapes of the kinetic energy spectra. Their values depend on the chosen 
fit range. If the fit range is increased it is very probable that, to obtain a good 
fit, also the multitude of temperatures has to be increased. This is the reason 
why the upper energy boundary in case of the vr" spectra was limited to 750 
MeV. The other reason is that the exact shape of the kinetic energy spectra at 
large kinetic energies, when the cross section has dropped by several orders of 
magnitude, depends on an accurate estimate of the background and on the way 
pions are selected from the data. The identification of pions with large energies, 
i.e. large momenta, requires to determine these momenta with high precision by 
the tracking algorithm. Above 750 MeV the standard deviation of the lab mo- 
mentum in forward direction becomes larger than 7% ^ and this is considered 
inadequate to identify such rare pions with the required precision. Nevertheless, 
the systematic variations of Ti^^ and T^^tt with system mass and energy allows 
one to make important conclusions. If the observed pions are assumed to be ex- 
clusively remnants from the decay of baryon resonances, then the increase of Ti ,r 
and Th^TT with increasing energy and mass is related to the change of the number 
of resonances involved and on their excitation profile at pion freeze-out times. In 
principle the pion energy spectra can be used to derive information about these 



quantities. This line of analysis shall be followed in a forthcoming paper p4[. On 
the other hand if the high-temperature component is interpreted as due to pions 
thermalized in the nuclear medium p6|, the values of Th^r allow one to deduce 



the fireball temperatures at pion freeze-out times. This line of analysis is followed 
in ref. However, one should not forget the problems which are caused by the 
quantum nature of pions and which were mentioned in subsection 3.2.1. 
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